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ABSTRACT 

Photoionization modeling of the low-ionization broad absorption lines of certain 
quasars, known as FeLoBALs, has recently revealed the number density of the wind 
absorbers and their distance from the central supermassive black hole. From these, the 
feedback efficiency of the quasars can in principle be derived. The implied properties 
of the FeLoBALs are, however, surprising, with the thickness of the absorbers relative 
to their distance from the black hole, AR/R, as small as ~ 10 . Such absorbers 
arc unlikely to survive the journey from the supermassive black hole to their inferred 
location. We show that the observed FeLoBAL properties are readily explained if they 
are formed in situ in radiative shocks produced when a quasar blast wave impacts a 
moderately dense interstellar clump along the line of sight. This physical picture dif- 
fers significantly from the thin shell approximation often assumed, and implies outflow 
rates, kinetic luminosities and momentum fluxes that differ correspondingly, in some 
cases at the order of magnitude level. Using the radiative shock model, we estimate 
the ratio of the outflow kinetic luminosity to bolometric luminosity for three bright 
FeLoBAL quasars in the literature. We find i^k/^boi ~ 2 — 5% (and corresponding 
momentum fluxes P s» 2 — 15 L^i/c), similar to what is adopted in models reproduc- 
ing the M — a relation. These outflow properties are also comparable to those recently 
inferred for molecular outflows in local ultra-luminous infrared galaxies, suggesting 
that active galactic nuclei are capable of driving such outflows. 

Key words: cosmology: theory - galaxies: active, evolution - quasars: general, ab- 
sorption lines 



> 1 INTRODUCTION 



X 



The broad absorption lines (BAL) seen in up to 40% of 
quasars (e.g., Trump et al. 2006 Dai et al. 2008), and 



blueshifted by up to ~ 0.1c relative to systemic, are clear sig- 
natures of outflows driven by active galactic nuclei (AGN). 
Such outflows are of particular interest in the context of the 
measured correlations between galactic-scale properties and 
those of the supermassive black holes (SMBH) they host, 
such as the M — a relation (e.g., |Ferrarese & Merritt||2000 



Gebhardt et al.||2000| |Tremaine et al.||2002[ ). Indeed, it has 



been argued that feedback by accreting SMBHs could regu- 
late their growth and produce these relationships (e.g., [Silk 



k Rccs 1998; Wyith e~fcLoeb||2003| [Pi Matteo et al.|2005[ ). 
While these arguments are broadly compelling, it remains an 
open question how the energy and/or momentum released 



by the SMBH couples to the surrounding galaxy, quenches 
black hole accretion, and possibly truncates star formation. 
Observational constraints on the feedback provided by AGN 
are therefore extremely valuable, but have been scarce so far. 

Progress on this front has recently been achieved 
through the detailed study of a subset of BAL systems, 
for which photoionization modeling has allowed inference 
of their column density Ah and radius R from the SMBH, 
in addition to their velocity v (e.g., Moe et aL"j|2009 Dunn 



et al. 2010 Bautista et al. 20101. Because these systems 



show absorption from low-ionization states of iron, specifi- 
cally Fell, they are known as FeLoBALs]^] Assuming a thin 
shell geometry, these measurements have been interpreted 
as providing estimates of the mass outflow rates, M out , and 
kinetic luminosities, of the quasar outflows. These stud- 
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1 Hel, Sill, and Nill ions also provide powerful diagnostics, when 
present. 
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Table 1. Sample of FeLoBALs and their Inferred Properties 



QSO 


z 


-^bol 1 


M BH 2 


log7V H 


logn e 


log[/ H 3 


T 


R 4 


AR 5 


V 






erg s — 1 


M 


cm 


cm - 3 




K 


kpc 


pc 


km s — 1 


SDSS J0838+2955 a 


2.043 


10 47.5 


3 x 10 9 


20.8 


3.8 


-1.9 


w 10 4 


3.3 


0.04 


4,900 


SDSS J0318-0600 b 


1.967 


10 47.7 


4 x 10 9 


20.1 


3.3 


-3.0 


ft 10 4 


5.5 


0.02 


4,200 


QSO 2359-1241° 


0.858 


10 46.7 


4 x 10 8 


20.6 


4.4 


-2.4 


ft 10 4 


1.3 


0.005 


1,400 



10 8 M (L bol /1.3 x 10 46 ergs" 1 ). 



1 Bolometric luminosity of the quasar. 
2 SMBH mass, assuming that the bolometric luminosity equals the Eddington luminosity, i.e. Mbh 

3 Ionization parameter Uh = Jn/oin, where <I>h is the hydrogen ionizing photon flux. 
4 Distance of the absorber from the central source. 
5 Thickness of the absorber, defined as AR = 1.2Nn/n e , assuming that the system is of solar composition and fully ionized. 
^Component c \Moe et al.|2009| . 'Component i ( |Dunn et al.|2010| . c Component e jKorista et al.|2008| [Sautista et al.|2010[ |Arav|2010[ ). 



ies of a handful of FeLoBALs indicate that in some cases the 
kinetic luminosity of the AGN outflow can exceed f = 1% 
of the bolometric luminosity (e.g., Moe et al. 20091. Such 
values are not far from the value / = 5% thought to be suf- 
ficient to explain the normalization of the M — a relation 
(e.g., Di Matteo et al. 2005; see also DeBhur et al. 2011). 

The physical properties of the FeLoBALs from which 
the kinetic luminosities and mass outflow rates are esti- 
mated (Table [T| are surprising and differ significantly from 
the interpretation of other BAL observations. In particular, 
whereas classical CIV BALs are inferred to be located at 



distances R < 1 pc from the SMBH (e.g., Capellupo et al 
2011[ ), the FeLoBALs in Table [l] have R > 1 kpc. Further- 
more, the small inferred FeLoBAL sizes AR ~ 0.01 pc im- 
ply large ratios R/ AR ~ 10 and, as we will show, lifetimes 
~ 10~ 2 of the time it would take such absorbers to travel 
from the SMBH to their inferred locations, tfl ow ~ R/v. To 
determine the reliability of these observational estimates of 
the efficiency of AGN feedback, it is important to have a 
more detailed physical picture of the FeLoBAL absorbers. 
Such a model is in particular necessary to determine how 
global parameters of the outflow can be derived from the 
absorption measurements. 

In this paper, we argue that the FeLoBALs must be 



formed in situ at large radii ((2.1 1. We show that the inter- 
action of a quasar blast wave with a moderately dense inter- 
stellar clump along the line of sight can result in radiative 
shocks, and that cool material entrained by the hot ambient 
flow can explain the properties of the observed FeLoBALs 



({2.21. This physical picture has important implications for 
the kinetic luminosities of the quasar outflows probed by 
low-ionization absorption lines (93J. We conclude in {|4]by 
discussing our results in the context of models of the co- 
evolution of black holes and galaxies. 

Before proceeding, a clarification of the terminology 
adopted here is in order. The different classes of broad 
quasar absorbers are usually defined in terms of their ob- 
servational spectroscopic properties. While these are indica- 
tive of their physical nature, there need not always be a 
one-to-one correspondence. In particular, only a handful 
of FeLoBAL absorbers to date have been studied in suf- 
ficient detail to extract their physical properties. In this 
work, we assume that the systems in Table [T] are represen- 



tative of FeLoBALs in general. Strictly speaking, it is con- 
ceivable that some FeLoBALs are physically distinct and 
formed differently than we propose (e.g., |Hall eFldTpnTj ). 
Furthermore, some members of the broader class of low- 
ionization broad absorbers (LoBALs), or even some high- 
ionization BALs (HiBALs), may form as we describe here. 
More detailed observations and photoionization analyses will 
be needed to determine exactly how widely applicable the 
model presented here is. 

Our model is based on the crushing of an interstellar 
cloud by the forward shock driven into the ambient medium 
by the supersonic quasar outflow. A review of the important 



cloud crushing physics can be found in Klein et al. ( 1994 1 



A key element is that the shocks driven into the over-dense 
cloud are much weaker than the forward shock in the am- 
bient gas, which allows them to cool and dust to survive 
({2.2|. A third shock involved is the reverse shock into the 



quasar wind itself. 
2 MODEL 

2.1 FeLoBALs are Created In Situ 

Quasar wind ejecta with velocity v take a flow time 

R « . I R 



til. 



ow ~ - w 3 x 10° yr 

v \,i kpc 



10, 000 km s" 



(1) 



to travel out to a radius R^ As the ejecta impact the in- 
terstellar medium (ISM) of the galaxy, the gas behind the 
shock is heated to a temperature 



(2) 



where fi is the mean molecular weight (assuming [i = 0.61 for 
a fully ionized gas of solar composition) and k is Boltzmann's 
constant. In this expression, we have assumed v = 3i) s h/4, as 
appropriate for the post-shock material for a strong shock. 



2 If the ejecta decelerate by a factor of a few along the way, the 
required time scale is correspondingly shorter, but this does not 



significantly affect our arguments. 
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n H pre , T f 



QSO blast wave encounters moderately 
dense ISM cloud. 




Shock wave propagates in cloud on 
crushing time tec, cloud is destroyed 
by K-H in tKH~20tcc, and is accelerated 

to ~Vsh in tdrag. 




At t>tKH, tdrag, original cloud is shredded 
into cloudlets traveling at ~vsh and 
compressed by hot post-shock gas. 



Figure 1. Schematic illustration of the formation of FeLoBALs in radiative shocks. 



For comparison, the sound crossing time of a cloud of size 
AT? is 



_ AR / AR \ ( T 



(3) 



so that it would have ample time to adjust its size. If the 
cloud has temperature Tbal, number density nj| AL , and is 
pressure confined by a medium at temperature T s h, that 
medium must have density 



hot Tbal bal 
1h = — n H 

ish 



(4) 



0.08 cm" 



/ Tbal 



VlO 4 cm" 3 J \10 4 KJ V 1-3 x 10 9 K 



Non-gravitationally bound cool clouds (such as those 
m Table [TJ moving in a hot medium are disrupted by the 
Kelvin-Helmholtz instability on a time scale 



„BAL \ 1/2 a D 

n w \ AR 



(5) 



where the velocity is relative to the hot medium, and the 
dimensionless factor can have values up to re w 10 when 
cooling is efficient (e.g., Mellema et al.|[2002 Cooper et al. 
|2009[ ), but not much more. For the fiducial parameters in 
equations |3| and Q, assuming that v = c s (T s h), equation 
|5| yields i K H ~ 630re yr. Since t K n <C tfi ow for the FeLoB- 
ALs in Table [l] these absorption systems would have been 
completely destroyed along the way if they had been pro- 
duced in the immediate vicinity of the SMBH, as is inferred 
to be the case for many classical, high-ionization BAL ab- 
sorbers. While the Kelvin-Helmholtz instability can in prin- 
ciple be suppressed by magnetic fields, the Rayleigh- Taylor 
instability generically cannot (Sto ne fc Gardiner|2007 1, and 
would destroy the clouds on a comparable time scale. Fur- 
thermore, thermal conduction from the hot medium would 
evaporate the cool clouds on a time scale 



^cvap 



6 x 10 3 yr 



AR 
0.01 pc 



7/6 



„BAL 

10 4 cm- 



n n 1 BAL 
10 s K cm" 3 



(6) 

-5/6 



(saturated heat flux regime; Cowie fc McKee]|1977 1 regard- 
less. 

Unlike the BALs directly accelerated near SMBHs, for 
example via disk winds (e.g., Murray et al. 1995 I, the FeLoB 



ALs in TableQ]from which mass outflow rates and kinetic lu- 
minosities have been estimated therefore very likely formed 
at the large radii where they are observed. 



2.2 Formation in Radiative Shocks 

At velocities > 1, 000 km s _1 , the post-shock gas cools pre- 
dominantly by free-free emission, losing its energy at con- 
stant pressure on a time scale 



10 s yr 



/ Tsh 



VlO 9 K 



1/2 



1 cm 



(7) 



where we assume that the shock is strong, so that the post- 
shock density is 4 times the pre-shock density n^ rc . A limit 
on the average number density swept up by the quasar out- 
flow can be derived using momentum conservation. 

Over a time scale ifl ow , the quasar injects a momentum 



p _ , Tjboi 

' — Jw ''flow 
C 



(8) 



into a wind. The factor / w can be larger than unity, for ex- 
ample in a confined shocked wind doing work on surrounding 
gas or if there are multiple photon scatters in an optically 
thick medium (e.g., Murray et al. 20051; alternatively, the 



mass inflow and outflow rates can in principle exceed the 
Eddington rate. If p pre (< R) is the average gas mass den- 
sity enclosed within a radius 7?, where the velocity of the 
outflow is v, 



AttR j 



p prc (< R)v. 



(9) 



Using n^ re (< R) = Xp pTC (< R)/m p (where X = 0.75 is the 
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hydrogen mass fraction), we obtain 
3-X/w Lhoi 



< c « R) 



4irR 2 cm p v 2 



0.02 cm"7. 



(10) 



10 47 erg s- 1 J \3 kpc 

2 



R 



5, 000 km s- 



At such densities, the post-shock gas would have no chance 
to cool in a flow time, making it impossible to form FeLoB- 
ALs at large distances via interaction with "average" ambi- 
ent conditions. 

Because FeLoBALs are rare (only ~ 1/1,000 of bright 
optical quasars show transitions suitable for photoionization 
analysis; e.g., Moe et aL]|2009 1, however, we may appeal to 
rare circumstances to produce them. Consider the situation 
illustrated in Figure [I] A moderately dense clump, for ex- 
ample analogous to a cold neutral medium cloud, lies along 
the line of sight. Let n^ 1 be the initial density of the cloud 
and aj its radius. The cloud is impacted by a blast wave 
from the quasar traveling at velocity v B ^ » 5,000 km s , 
as above. As the blast wave overtakes the cloud, a shock 



wave is driven into it at a velocity u s h, c ~ «shy "h"/^'' se ^ 
by the requirement that the post-shock gas in the cloud be 
in pressure equilibrium with the surrounding medium (e.g., 
McKee fc Cowie|1975|. The corresponding post-shock tem- 



perature is then T s h, c = T B h(ri%[ e /n^ 1 ), and we define a cloud 
crushing time as 



Cti 
Vsh,, 



I'sh 



(11) 



Identifying the cloud size A_R = 2ai and v = 3« s h/4, equa- 
tion (JsJ implies that the cloud is destroyed in £kh = 8fti cc /3. 

At the same time, ram pressure from the hot exter- 
nal medium accelerates the cloud, which becomes co-moving 
with the hot shocked gas on a drag time 



tdt 



4 ai n. 



3 v 4n% e ' 



(12) 



i.e. the time necessary for the hot flow to sweep the cloud 
with a mass equal to its own. The requirement that the cloud 
be accelerated before it is torn apart is 



tdi 



< 1 «■ 



(13) 



We further require that the cooling time behind the shock 
driven into the cloud be less than a cloud crushing time, so 
that the shock is radiative and leaves cold material behind. 
Such radiative cooling also prolongs the cloud survival with 
respect to the Kelvin-Helmholtz instability (e.g. 
et al.|2002 Cooper et al.|2009 \ , allowing values k 



Mellema 



10. For 

10 s < T < 4 x 10' K, the isobaric cooling time for ionized 
gas is 



fcooi ~ 1.5 x 10 yr 



sh,c 



10 6 K 



1.6 



4"h 
10 cm- 



(14) 



where we have adopted the approximation to the cooling 
function from |McKee fc Cowie] ( |1977[ ) , and the formula ac- 
counts for a factor of 4 compression through the strong 



shock. Combining the expressions above, 



^cool 
tec 



■I'sh 



5,000 km s" 1 

2.1 

18 cm" 3 



0.02 enr 



(15) 



1.4 x 10 20 cm" 2 



where N^ 1 = 2ain^' is the column density of the original 
cloud. In this formula, the fiducial values for n^ re and n^ 1 
were chosen to satisfy fdrag ~ tra with k = 10 (c.f., eq. 
| X 3 P - i.e. to ensure that the cloud is accelerated before being 
shredded. It follows that column densities TV^j' 1 > 10 20 cm -2 
are necessary for the post-shock gas with v s h ~ 5, 000 km 
s _1 to cool and produce conditions suitable for FeLoBALs. 

For t > tdrag, £kh, the main cloud has been shredded 
into multiple cloudlets traveling at v ~ ^ s h- Since t coo i < t cc , 
these cloudlets cannot maintain pressure balance with the 
surrounding medium and are therefore compressed. These 
photoionized cloudlets cool to T as 10 4 K, where the cooling 
function is strongly suppressed, and reach densities up to 
(in pressure equilibrium with the hot gas) 



4nT 



10 4 K 



(16) 



3,000 cm" 



0.02 cm- 



Vsh 

5,000 km s" 1 



The predicted densities, column densities]^] and tempera- 
tures of these cloudlets are in excellent agreement with the 
inferred properties of the FeLoBALs summarized in Table 

in 

Furthermore, this picture naturally explains the obser- 
vations of multiple absorption components along the line of 
sight with approximately the same distance from the SMBH 
(e.g., Bautista et al.|[2010 |, as some cloudlets may only be 
partially accelerated. The supra-thermal line widths mea- 
sured (hundreds of km s~ ) could also result from a combi- 
nation of blended components and velocity shear during the 
shredding process. 

In our model, the fact that FeLoBALs arise along par- 
ticular sight lines with over-dense material is also in agree- 
ment with their host quasars being redder than average 



(even among BAL quasars; e.g., Sprayberry & Foltz 1992 
|Becker et aLl|2000| |Urrutia et al'||2008[ ). It is also possible 
that dust is present in other clouds along the line of sight, 
which in general could be seen in absorption at different ve- 
locity offsets. However, since strong absorption at the galaxy 
systemic redshift is often not seen in FeLoBALs, this sug- 
gests that the dust responsible for the reddening is part of 
the FeLoBAL gas, and therefore that the dust survived the 
cloud crushing shocks. We note that the intra-cloud shocks 
have velocities much lower than the quasar blast wave ow- 
ing to the cloud over-density. Even for a quasar blast wave 
with velocity v s h = 5, 000 km s" 1 , for the maximum density 



3 In the radiative case, the original cloud does not undergo signif- 
icant lateral expansion (as it would adiabatically; e.g. Kle in et al.| 
|1994[ |. At least before they mix significantly with the hot medium, 
the integrated column density of the cloudlets is therefore sim- 
ilar to that of the original cloud, provided that their collective 
covering factor is preserved during the shredding. We therefore 
identify N^ 1 



N* AL . 
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contrast allowed by equation ( 13 1 with k = 10, the cloud 
crushing shocks have velocities w s h,c ~ 170 km s _1 . At such 
velocities, it is reasonable for at least some dust to survive 



in the FeLoBAL gas (e.g., |McKee||1989| >. Indeed, although 
the dust is likely to be eventually destroyed with the cloud 



itself, the short lifetimes of order a few 1,000 yr ({2.1 \ im- 
ply that there is little time for dust sputtering processes to 
operate while the cloud is present. 

2.3 Selection Effects and Predicted Trends 

In reality, the hot flow likely entrains a wide spectrum of 
clouds. In the scenario outlined above, selection effects ex- 
plain the relatively narrow range of observed FeLoBAL prop- 
erties. 

Absorption by low-ionization species such as Fell, Hel, 
and Sill useful for photoionization modeling selects ioniza- 



tion parameters Uh 



lo- 



ur 



as well as tempera- 



tures sufficiently low to avoid collisional ionization of those 



species. As equations ( 13 \ and ( 15 \ show, cool gas remains 



behind the blast wave and is accelerated to a velocity com- 
parable to the hot flow only for clump column densities 



AC 



> 10 cm 



"sh 



5, 000 km s- 



(17) 



Interestingly, there is evidence that FeLoBALs located closer 
to SMBHs (~ 1 — 10 pc vs. > 1 kpc) have column densities 
higher by factors ~ 10—100 (e.g., |Dunn et al.|2010| ), a trend 
in agreement with our model if the blast wave decelerates 
with radius 

Our model also naturally explains why the FeLoBAL 
absorbers for the bright quasars in Table [I] are all located 
at relatively large radii R ~ 1 — 3 kpc. Indeed, the obser- 
vational selection of a narrow range of ionization param- 
eter !7h oc Lboi/ R 2 nH defines a relationship between the 
quasar luminosity Lboi and the FeLoBAL radius R. Since 
the FeLoBAL gas density nu is set by pressure equilibrium 



as in equation ( 16 1, this relationship in fact implies that the 



FeLoBAL absorbers in luminous quasars would be typically 
found at large radii. 



3 INFERRING THE FEEDBACK EFFICIENCY 

The physical picture of FeLoBALs developed here allows us 
to relate their observed properties to those of the underlying 
quasar outflows. 

Consider first the common assumption that the outflow- 
ing material is confined to a partial, cold, thin shell that has 
a covering factor Q and a column de nsity equal to that of 
the FeLoBAL, Af| AL (e.g., |Arav|2010[ ). Since the mass of the 
shell Mshoii = 4:TvR 2 N^ Alj mpQ/X, the mass outflow rate is 
given by 



Af 8 , 



8ivRNi Ah m p v 
X 



o. 



(18) 



4 Since not all cloudlets need be fully accelerated to the hot flow 
velocity, their velocity may not always reflect the blast wave ve- 
locity. Furthermore, the criterion in equation \17\ is strengthened 
for values of n^ 1 /n^ c smaller than fiducially assumed in equation 
| |15[ l . There is therefore not necessarily a well-defined relationship 
between Nh and R, or TVh and v. 



where v = R. The corresponding kinetic luminosity and mo- 
mentum flux are then 



and 



MshellW- 



(19) 



(20) 



Motivated by the observed incidence of general (mostly 
high ionization) broad absorption lines in quasars, these ex- 
pressions have usually been applied with a covering factor 
n ~ 0.2 (e.g., |Moe et al.|2009||Dunn et al.|20lo| |. We return 
to this point below, where we argue that this may generally 
not be a valid assumption for FeLoBALs. 

In our picture, FeLoBALs are rare tracers of an under- 
lying shock-heated medium, fortuitously encountered along 
the line of sight. A more appropriate estimate of the mechan- 
ical properties of the outflow can be obtained by replacing 
the column density of the FeLoBAL by that of the hot gas 



in equation (18 1, and similarly for the velocity and covering 
factor, i.e. 



8tt RN& ot m p v hot 
X 



(21) 



The shell approximation remains reasonable, since at large 
radii most of the hot gas is swept up, shocked ambient ISM. 

The properties of the hot flow required in equation |21| 
are not directly measured, but can be inferred by making 
assumptions motivated by our model. One subtlety is that 
in a realistic outflow, the ambient medium is separated from 
the direct quasar wind by a forward shock, a contact discon- 
tinuity, as well as a reverse shock. We discuss this in detail in 
the Appendix but it does not change the conclusions drawn 
here using simpler arguments. 

We begin by assuming that the FeLoBAL gas is fully 
accelerated by the hot gas, v « Vhot, and in pressure equi- 
librium with it: 



hot 

n H 



n BAL /10 4 K 

V 2hot 



(22) 



As before, T bot » T sb (v sb = 4v/3), where T sh is given by 
equation Q. The column density of the hot gas can be es- 
timated as 



Nk ot w n%°(< R)R- 



(23) 



If the FeLoBAL is located in the shocked ambient medium 
and close to the forward shock (see the Appendix), then 
riH 0t ~ 4nj^ rc . Under these assumptions, 



1 W BAL ( R \ f lO* K 

4 Nh I AT? J Vt^ 



(24) 



where AR is the thickness of the FeLoBAL absorber (as be- 
fore). Written this way, all the quantities necessary to evalu- 
ate equation (21 1 are provided by photoionization modeling 
of the FeLoBAL (e.g., Tablejl]), except That which is derived 
from v, and Qhot- The energy and momentum fluxes of the 
outflow follow straightforwardly given Afhot and Vh„t. 



21 



is not 



It is important to stress that Mhot in equation 
the intrinsic mass outflow rate in the quasar wind. Instead, 
Mhot is the rate at which the ambient ISM material is being 
swept up by the quasar outflow, i.e., it is the galaxy-scale 
outflow rate of shocked gas. In the model of FeLoBALs ad- 
vocated here, the observations cannot directly constrain the 
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Table 2. Galaxy-Scale Mass Outflow Rates, Kinetic Luminosities, and Momentum Fluxes Inferred from the FeLoBALs in Table [l] 



QSO 


M shell (O = 0.2) 


^shcll 


^jshcll 


Mhot(^hot — 1) 


^hot 




pQSO 




M yr- 1 


%Lbol 


Lbo\/c 


M Q yr- 1 


%-Lbol 




%L bol 


SDSS J0838+2955 


560 


1.2 


1.5 


1,000 


2.2 


2.7 


3.5 


SDSS J0318-0600 


150 


0.2 


0.2 


1,100 


1.2 


1.8 


2.0 


QSO 2359-1241 


36 


0.05 


0.2 


2,400 


3.1 


13.4 


5.0 



Note. — Values assuming that the outflow is confined to a cold thin shell (eq. |18| are labeled 'shell.' Values assuming that the 
FeLoBALs arise in radiative shocks (eq. \21\ and are entrained by an energetically-dominant hot flow are indicated by the label 'hot.' For 
these, the FeLoBAL gas is assumed to be co-moving and in pressure equilibrium with the shocked ambient medium. A covering factor 
Q = 0.2 is assumed for the cold thin shells for comparison with previous work, but a more physically plausible covering factor Qhot = 1 is 
adopted for the hot gas at large radii (q3f . The last column, E^ SO , lists the kinetic luminosities of the small-scale quasar wind inferred 
assuming a self-similar adiabatic wind solution (see the Appendix). If the FeLoBAL gas is not fully accelerated to the velocity of the hot 
gas, Mhot should be interpreted as an upper limit, and E^ ot and E^ SO as lower limits. E^ SO provides a more complete estimate of the 
quasar wind energetics than E^ ot , which neglects the fraction of the energy thermalized in the hot flow. 



intrinsic quasar mass outflow rate because most of the mass 
that is pressure confining the BAL gas is likely swept-up 
ISM material, rather than the direct quasar wind. However, 
the observations do directly constrain the galaxy-scale out- 
flow rate Mhot and the corresponding energy and momentum 
fluxes ££ ot and P hot . Note that E^ ot corresponds only to the 
bulk kinetic energy of the hot gas, neglecting a fraction ~ 1 
that is instead thermalized. The derivation of the small-scale 
quasar wind kinetic luminosity £^ so for a self-consistent 
outflow solution in the Appendix, on the other hand, ac- 
counts for all the energy. Accordingly, E\ 



qso 



2E? 



In Table [2] we give the galaxy-scale mass outflow rates, 
kinetic luminosities, and momentum fluxes derived from the 
FeLoBAL systems in Table [T] evaluated using both the cold 
thin shell approximation (eq. |18[ ), and assuming that the 
FeLoBAL is only a tracer of the underlying shocked medium 



as above (eq. 21 1. For the former case, we fiducially adopt 
Q = 0.2 for comparison with previous work. For the lat- 
ter, we instead fiducially assume Qhot = 1. A covering fac- 
tor Qhot ~ 1 seems more physically plausible at kilo-parsec 
radii, because in the absence of significant obstacles the hot 
bubble will tend to become more spherical as it expands. 
This applies even if the energy is, for example, injected in 
the form of high-ionization BAL systems with ~ 0.2 from 
an accretion disc wind at R < 1 pc. Momentum conservation 
actually suggests that the outflows traced by the FeLoBALs 
in Table [T] are propagating along sight lines of low gas den- 
sity, consistent with a spherically symmetric geometry and 
possibly owing to an earlier phase of feedback which cleared 
most of the gas out. A dense galactic disc might impede the 
outflow along certain directions, but should have a minor 
effect provided it is thin. We however cannot exclude the 
possibility of viewing the quasars along clear sight lines by 
chance. Thus, there is significant uncertainty in the normal- 
ization of the energetics in Table[2]due to the covering factor 
assumption. 



normalizing to the same covering factor]^] Indeed, this must 
generically be the case for fully accelerated FeLoBALs, since 
then 



N hot 



"hot 



(25) 



and full acceleration requires 

Nu" > Ni Ah (i.e., that the 
cool cloud be swept with a mass equal to its own). Another 
bound on the ratio of the column densities of the cold and 
hot gas can be derived in the context of our model by com- 
bining the momentum conservation and cooling conditions 



(eq. 10 and 17 1 



jyhot 



Lb 



/ w V 1047 erg s -1 



R 



3 kpc 



(26) 



5, 000 km s- 



This expression is, again, in good agreement with the re- 
sults in Table [2] Indeed, for SDSS J0838+2955 (v = 4,900 
km s _1 ) Afshdi ~ Mhot for the same covering factor, as com- 



bining equations (|25p and ( 26 1 indicates should be the case 



However, the steep dependence of equation ( 26 1 on v implies 
that M 6heU for QSO 2359-1241 (v = 1,400 km s" 1 ) can be 
much less than Mhot, and in fact it is found to be. 

Although we have so far assumed that the FeLoBAL 
gas is fully accelerated, i.e. co-moving with the hot flow, this 
may not always hold. Provided only that the cool absorbers 
are pressure confined, Mhot oc v^ t , E^ ot oc Uhot, and P hot 
is independent of fhot- In particular, a useful lower limit on 
the kinetic luminosity of the outflow is obtained (for a given 
f2hot) if the FeLoBAL velocity v ^ Mhot, while the derived 
Mhot is then a lower limit. The values in Table [2] should 
more generally be interpreted this way. 

FeLoBALs with properties identical to those in our 
blast wave model can also be produced in a smooth, su- 
personic wind (with a density profile oc R~ 2 and velocity 



As Table [2] shows, the mass outflow rates inferred for the 
hot flows (and consequently E^ and P) are systematically 
higher than those implied for the cold thin shells, even after 



5 This is violated by a factor < 3 for SDSS J0838+2955, but this 
is only marginally significant at the level of our estimates. It is 
also possible that the absorber is not fully accelerated. 
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Win) entraining an interstellar cloud, since the cloud crush- 
ing then proceeds similarly (e.g., Klein et al.|1994 l. In this 
case, the cold gas reaches a terminal velocity v satisfying 



{v/viu) < 1, where N£ is the column density 



through the wind at radius R. This "pure wind" scenario is 
less physically plausible since FeLoBALs have velocity off- 
sets ~ 1,000 — 5,000 km s _1 , compared to accretion disc 
wind velocities ~ 0.1c = 30,000 km s _1 . This suggests that 
the quasar wind has decelerated significantly along the line 
of sight by sweeping up ambient ISM. 



4 DISCUSSION 

We have argued that the FeLoBAL absorption systems seen 
in a small fraction of quasars are physically distinct from the 
larger class of broad line quasar absorbers. The majority of 
high-ionization BALs are located at radii R < 1 pc from 
the accreting supermassive black hole and are likely directly 
accelerated by it, possibly in an accretion disc wind (e.g., 
Murray et al.|1995 1. The kinematics of most LoBAL quasars 
also appear to differ from those of FeLoBALs: |Voit et al.| 
( 1993 1 stress that in the LoBALs they study (which do not 
show Fell absorption) the Mg II and Al III features are 
narrower (by factors of several in some cases) than the C IV 
troughs, and the Mg II and Al III features tend to lie at the 
low-velocity ends of the C IV trough. This contrasts with 



the FeLoBALs discussed by Moe et al. (20091, Dunn et al 



(20101 and Bautista et al. (20101, where the C IV troughs 



have similar total velocity width to that of Mg II and Al III, 
as well as similar velocity structure, in which the individual 
features of different ions line up in velocity space. The latter 
behavior strongly suggests that the C IV absorbing gas is 
co-spatial with the Mg II and Al III gas in FeLoBALs, unlike 



the C IV gas in the Voit et al. ( 1993 1 objects 



Consistent with our argument that the FeLoBAL sys- 
tems are physically distinct from other BAL systems, we 
have shown that the FeLoBALs form in situ in the ISM of 
the host galaxy. Indeed, when a quasar blast wave impacts 
a sufficiently dense clump along the line of sight, the cloud 
crushing shocks are radiative and leave cool material behind. 
We derived the conditions under which (1) the post-shock 
gas cools and (2) the cool gas is accelerated to a velocity 
> 1,000 km s -1 : the required properties are in excellent 
agreement with those inferred for FeLoBALs. 

This new physical picture of FeLoBALs has important 
implications for measurements of the efficiency with which 
quasars convert their luminosity into mechanical feedback 
on galactic scales. Unlike the cold thin shell approximation 
often assumed (e.g., Arav 20101, in our picture FeLoBALs 



are only rare tracers of an underlying outflow. In particular, 
the observed FeLoBAL properties constrain the properties of 
the ambient ISM shock heated by the quasar blast wave. We 
find that the mass outflow rate, kinetic luminosity, and mo- 
mentum flux of the hot flow can be larger, sometimes by an 
order of magnitude or more, than those previously estimated 
for the cold gas only. In Table [2] we give the values derived 
for three bright quasars in the literature, using both the cold 
thin shell approximation and in the context of our radiative 



shock model. In our model, E® so « 2 
responding momentum fluxes P h 
case of QSO 2359-1241, the kinetic luminosity E\ 



5% Lboi, with cor- 
— 15Lboi/c. In the 



Lboi is ~ 100 x larger than implied by the cold thin shell 
approximation. 

Our inferred values for the kinetic power in quasar out- 
flows on galactic scales are comparable to those assumed in 
models successful in reproducing the M — a relation (e.g., 
Wyithe fc Loeb|2003| [Pi Matteoet al.|2005[ ). Moreover, for 
a given covering factor of hot gas fihot ~ 1, our inference 
of E^ 80 is a lower limit if the cold BAL gas is not fully 
accelerated to the velocity of the hot gas. The momentum 
fluxes we find for FeLoBALs in Table [2] are also compara- 
ble to the values P ~ 5 — 10Lboi/c that appear sufficient 
to drive powerful galaxy-wide outflows in simulations that 
model black hole feedback by injecting quasar winds into 
the surrounding ISM ( |DeBuhr et al.|2011 1. 

The mass outflow rates and energy and momentum 
fluxes inferred in Table [l] apply on the galactic scales where 
the absorbers are detected (R ~ kpc). These need not corre- 
spond to the small-scale mechanical properties of the quasar 
wind. In particular, galactic-scale mass outflow rates are in 
general higher than those of the direct quasar wind because 
of swept up material: the large values in Table [T] therefore 
need not imply super-Eddington accretion disc winds. 

The FeLoBALs in the three bright quasars studied here 
may trace a particular stage of feedback. In fact, the obser- 
vations of cool absorbers moving at v ~ 5,000 km s" 1 at 
R ~ 3 kpc from the black hole imply relatively clear sight 
lines. Otherwise, the quasar wind would have been deceler- 
ated more significantly by sweeping up ambient ISM. One 
possibility is that we are witnessing late accretion events, 
with the outflows propagating into a tenuous ISM that has 
already been mostly evacuated in an earlier phase. The pow- 
erful molecular outflows recently detected in local ultra- 
luminous infrared galaxies (ULIRGs) may provide evidence 
of such gas clearing (e.g., Feruglio et al.|2010 Fischer et al. 
2010 Sturm et al.||20lT l, though star formation may also 
contribute to driving these winds (Chung et al. 20111. In- 



triguingly, the kinetic luminosities and momentum fluxes of 
the ULIRG molecular outflows are similar to those we infer 
for the FeLoBAL quasar outflows, suggesting that AGN are 
indeed capable of driving the ULIRG outflows. 

Independent constraints on the evolutionary stage of 
FeLoBAL quasars, for instance from their broad band spec- 
tral energy distributions, would be valuable in distinguishing 
the late feedback possibility from that of a chance geometri- 
cal effect. Subject to the significant limitations of small and 
heterogeneous samples, recent observations appear to sup- 
port the hypothesis that at least some FeLoBAL quasars are 
observed toward the end of their starburst and in the last 



stages of blowing out their surrounding ISM ( Urrutia et al 
[2^fj|{[Farrah et al.|2007|[20Tot . Since much of the black hole 
growth is predicted to occur while the AGN is obscured, 
appearing in the optical only after it has blown out the ob- 
scuring material (e.g., [Hopkins et al. 2005), this late stage 
of feedback picture is consistent with the "youthful" quasar 
hypothesis for FeLoBALs ( THall et al.|2002 i. 



QSO 



5% 
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APPENDIX A: RELATIONSHIP BETWEEN 
THE SHOCKED AMBIENT MEDIUM AND 
THE QUASAR WIND 

In ij3] we inferred the mechanical properties of the quasar 
outflow based on those of the hot, shocked ambient medium 



at R. In realistic outflow solutions (see Fig. Al |, the quasar 



wind is separated from that medium by a reverse shock at 
R BW and a contact discontinuity at R c . We denote the radius 
of the forward shock by R s . 

An important question is how E^ ot relates to the ki- 
netic luminosity of the wind in the immediate vicinity of the 
quasar, E® SO . For some purposes, such as comparing with 



Di 



galaxy simulations that inject energy near the BH (e.g 
Matteo et al. 2005), this quantity is more directly relevant 
than the hot flow on galactic scales. From FeLoBALs at kilo- 
parsec radii, this can only be done assuming a global outflow 
solution. 

As we showed (eq 



10 1, both the reverse and forward 



shocks are likely to be adiabatic for the FeLoBALs in Ta- 
ble [l] For the case of a quasar wind with constant E® so = 
^Mi n vf n (where the subscript 'in' denotes the initial proper- 
ties of the wind) propagating in an uniform ambient medium 
of density p pre , Koo & McKee (19921 showed that 



Rs = 0.; 



\ Pprc J 



3/5 



(Al) 



in this regime. This solution incorporates the shock jump 
conditions and the pressure balance requirement at the con- 
tact discontinuity in the dimensionless pre-factor. It also as- 
sumes that the outflow has reached radii sufficiently large 
that it is self-similar and no longer depends on the details 
of the energy injection. 

Taking the derivative, time can be eliminated from 



equation ( Al | and it is straightforward to show that 



■ 7 QSO 



8j 



r>2 p3 



(A2) 



Since the FeLoBALs are formed upon impact with the for- 
ward shock in our model, and i C va P <§C £flow for the sys- 
tems in Table [T] we can usually assume R ~ R E . In general, 
v 3.R s /4 and if the absorber is fully accelerated, then 
v w 3i? s /4. To estimate the density of the pre-shock ambi- 
ent medium, it is necessary to assume that the FeLoBAL gas 
is in pressure equilibrium with the gas behind the forward 
shock: 



ft* 



(A3) 



10 4 K 

Under these assumptions, which are essentially the same 
that were necessary to estimate the properties of the hot flow 
in 351 



E? SO < 2 x 10 46 erg s" 1 



BAL 

10 4 cm- 3 



R 



3 kpc 



(A4) 



s 5,000 km s- 

with equality holding in the case of full acceleration. This 



n H pre 



ambient " 
* medium^ 



QSO^ 




shocked 
wind 



sw 




Figure Al. Quasar outflow structure. The shocked ambient 
medium (bounded by _R S ), in which the FeLoBALs are envisioned 
to form, is separated from the direct quasar wind by a reverse 
shock at i? aw and a contact discontinuity at R c . 



can be compared to the kinetic luminosity of the hot flow 



estimated using equation (21) 



EC 



' r QSO 



0.6fi h 



(A5) 



a constant ~ 1. This near equality is not surprising since we 
assume that the outflow is adiabatic, so that the energy in 
the flow is preserved with radius. It does however formally 
confirm that the kinetic energy of the quasar wind can be 
estimated to a good approximation based on the proper- 
ties of the hot flow probed by FeLoBAL observations. An 
order unity difference between E^ ot and Els 80 is expected, 
as the derivation of E^ ot ignores the fraction of the quasar 
wind kinetic luminosity that is converted into thermal en- 
ergy, rather than bulk motion (§31. 



Since the self-similar solution in equation ( Al I depends 
only on E? so and Ppi- C , it is however not possible to sepa- 
rately infer M[ n or V[ n from the large-scale properties of the 
outflow. 
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